Regular and Chaotic Quantum Dynamic in Atom-Diatom 

Reactive Collisions 

Ashot S. GevorkyaEE] 
IIAP/IAPP NAS of Armenia, P. Sevak St. 1, Yerevan 375014, 

Alexander V. BogdancrvQ 
Institute for High Performance Computing and 
^ - Information Systems, St. Petersburg, Russian 

O ' Gunnar Nymarjj] 

> 

A new micro-irreversible 3-D theory of quantum multichannel scattering in the 
three-body system is developed. The quantum approach is constructed on the gen- 
1"**- ! erating trajectory tubes which allow taking into account influence of classical non- 

integrability of the dynamical quantum system. When the volume of classical chaos 
CNj ' in phase space is larger than the quantum cell in the corresponding quantum system, 

quantum chaos is generated. The probability of quantum transitions is constructed 
fS| . for this case. The collinear collision of the Li + (FH) — > (LiF) + H system is used 

for numerical illustration of a system generating quantum (wave) chaos. 
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INTRODUCTION 



In the early stage of quantum mechanics development A. Einstein asked a ques- 
tion that have attracted close attention several decades later [l[ . The question was: 
what would be the analogue of a classical chaotic system in quantum mechanics? In 
particular he pointed to the three-body system, which in general is well known to 
have a chaotic nature. 

In an effort to formulate and obtain the solution of the problem of quantum 
chaos, M. Gutzwiller tentatively divided all the existing knowledge of the dynamics 
of physical systems into three areas 0] : 

1. Regular classical mechanics (R area); 

2. Classical chaotic system or dynamical Poincare system (P area); 

3. Regular quantum mechanics (Q area ) . 
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The mentioned connected by certain conditions. Thus, Bohr's correspon- 

dence principle connects the R and Q areas, transferring quantum mechanics into 
classical Newtonian mechanics in the limit h — > 0. Areas R and P are connected by 
the Kolmogorov-Arnold-Moser (KAM) theorem. 

The general principle which can connect P and Q areas is not determined yet. Re- 
lated to the fourth, conditionally named the quantum chaos area Q c h, M. Gutzwiller 
mentioned that the "quantum chaos" conception is rather a puzzle than a well for- 
mulated problem. It is evident that the task formulated correctly in Q c h area is 
the most general one and under specific conditions must be transformed into the 
aforementioned limiting areas. 

Observation of chaotic phenomena in the spectroscopy of atomic nuclei (H, atoms 
0], molecules (5| and in billiard systems has stimulated a considerable interest 

in the quantum chaos problem in recent years. Irregular behavior of the wavefunction 
has been found in numerical calculations of quantum mechanical stadium billiard 
problem 0. It has been shown that the so called scars which were observed have 



classical trajectory characteristics [lOj. It has been known for a long time that 



classical models of chemical reactions exhibit chaos [Uj. It was shown that the 
mixing properties of chaotic dynamics observed in unimolecular reactions can be 
explained by some statistical laws [13]. Recall that one major motivation for the 
continued classical investigation of the reactive scattering problem is several 

kinds of experiments on waves, which have demonstrated the validity of the ideas of 



quantum chaotic scattering [15|-[17j. Atomic systems are quantum objects should 
thus be treated considering their quantum properties. 

The development of different semiclassical and mixed quantum-classical methods 
(see for example the detailed report (HI) can be considered as a natural extension of 
the classical trajectory study. This development has been motivated by the fact that 
the standard quantum approach is too demanding even for most few-body systems. 
For many problems various quasi-classical methods can give satisfactory results. The 
semiclassical methods, however, are restricted to relatively small systems. 

The problem of quantum chaos and its connection with classical nonintegrability 
was originally studied by the authors in the framework of a collinear three-body 
collision model [l9|. In the current article this approach is generalized to the 3D 
case. 



II. FORMULATION OF SCATTERING PROBLEM 

We will be interested in the three-body reactive scattering process A+ (BC) n — > 
(ABC)* — > (AB) m + C, where A, B, and C are atoms, n and m characterize the set 
of quantum numbers of diatomic states corresponding to initial (in) and final (out) 
scattering arrangements and (ABC)* denotes the resonance complex. Moreover 
m-A, n^B and mc are the masses of the particles and ta, tb and rc the column vectors 
describing their positions relative to an origin fixed in the laboratory system. The 
reactant arrangement is best described by mass scaled reactant Jacobi co-ordinates, 
while the product arrangement is best described by mass scaled product Jacobi 
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coordinates. For the reactant arrangement we can write 20|, |2jJ : 

A R,„, Qi n = A 1 r„, 



qo c 



(1) 



where R a and r a Jacobi coordinates of reactant (in) channel, moreover: A = 
[mi(l-m J 4/M)//i] 1 '' 2 , /i= [m^mBmc/M] 1 ' 2 , M = + m B + m C - In term 
of coordinates (qo a > Qia) the Hamiltonian of three-body system takes: 

H(q;P q ) = (l/2/i)Pj + y(g o ,gi,0), q = (qo, qi) = {q k }, k = 0,..,5. (2) 

Note that here and in the following we omit the channel index for simplicity. In (j2J) 
fj, and P q are the effective mass and moment of body system, (q = |q |, qi = |qi|, 9) 
characterizes the intrinsic coordinates, 9 is the angle between vectors qo and q x . The 
remaining coordinates (qs, (74,(75) are expresses via Euler angles. The interaction po- 
tential between all atoms V(q^,qi,9) depends on intrinsic coordinates. Recall, that 



a) 





'/i 



9i 



FIG. 1: a) Intrinsic Jacobi and local coordinate systems. The angle i9 is defined from 



cot # = b = y uiAmc I 'tubM . b) The reaction path is passing through the minimums 
of potential energy while the reaction coordinate can be an arbitrary smooth curve 
connecting (in) and (out) asymptotic channels. The lower shaded area is self-crossing 
region for natural collision coordinate (NCC) system associated to the reaction path curve 
while the upper one self-crossing region for NCC system is associated to the curve 



the coordinate systems needed for reactants and products are different [18j. This 
fact creates certain mathematical and computational complexities for the investiga- 
tion of multichannel scattering problem. The way to overcome it is to turn to special 
type of curvilinear coordinates, which are natural and suitable for description of two 
(or more) asymptotical states (in) and (out) simultaneously. For satisfying of this 
conditions in the collinear collision case was introduced smooth curve (coor- 
dinate reaction) which connected (in) and (out) asymptotical channels and along 
which was defined local orthogonal coordinates system (u, v) (see [22I. [23|). 
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In 3D case too we can introduce the curve Qw, along which NCC system is 



defined. In this case is defined on the plane (go, gi, # = 0) by expression [24 



(& = a/ {Qi ~ (leg) + Hi + qtqi Qe q < <f\ < 00, (3) 

where a and b are constants. In Eq. (j3J) q~ q and g+ are mass-scaled equilibrium bond 
lengths of molecules in the (in) and (out) channels respectively, a is an arbitrary 
constant, which is usually chosen to make the curve pass close to the saddle point 
of the reaction. The superscript c over go and <?i underlines the fact that the point 
(go, qf) lies on the curve. The limit (gi = 0, go — > oo) corresponds to the (in) state, 
while the limit (go = bqi = q\ cot d) corresponds to the (out) state. The movement 
along the curve Q^y is described by the coordinate u: 



u 



u - a/qf + bqi, q\ = q\- q eq , (4) 



where Uq is some initial point on the curve S$if. The inverse transformations between 
the two pairs of coordinates (go,gi) (u,v) are: 

q (u,v) = qn(u) —v sin 4>(u), qi(u, v) = q{(u) +v cos(j)(u), (5) 

where v is the distance from the curve Q^/. In Eq. (jSJ) the angle 4>(u) is determined 
by requiring orthogonality of coordinate system (u, v): 

dq^/dql = cot(j)(u), lim cot <f>(u) — cot i?. (6) 

(u,v=0) u^+oo 

Let us introduce the system of orthogonal local coordinates x = x(x°, x 1 , ...x 5 ) along 
the curve 0^/ using the transformations: 



x 



= u, x 1 = v, x 2 = f(u,v,9), x 3 = d Ui, x 4 = d u 2 , x 5 = d u 3 , (7) 



where function f(u,v,9) = \J (g ) 2 — 26g gi cos6 + b 2 (qi) 2 is mass-scale distance 
between A and B particles. In some part of the 3D Cartesian configuration space 
these equations determine a biunivocal mapping between the two intrinsic coordinate 
systems: {qo,qi,9} and x 1 , x 2 }. The set of coordinates (u±, 0^2,^3) are three 
Euler angles, which orient the three-body system in the space-fixed frame [25| , <i is 
some space-dimensional constant. 



A. Classical dynamics of three-body scattering system 

For the investigation of ergodic properties of conservative dynamical system the 
geodesic axes distribution method on Riemann surfaces had been originally applied 
in 26[ . Later this method has been used and developed in the investigations of 
the foundations of statistical physics (13]. The study of geodesic flow behavior 
on Lagrange surfaces provides an opportunity to observe important properties of 
classical dynamics systems [28| . 
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Consider the QD three-body classical problem on the Lagrange surface Sp: 

S P = {x; P 2 (u,v,9) = 2fi[E-U(u,v,0)] > 0}, (8) 

where E is the total energy and U(u,v,9) is the interaction potential of the three- 
body system. The metric on the surface Sp is introduced in conform-Euclidian 
form: 

(ds) 2 = ^gijdx'dx 1 , gij(u,v,9) = P 2 (u,v,9)Sij, j = 0,1,..., 5. (9) 

Now we can write the geodesic trajectory problem for the reduced mass fi: 

■i-L ■ l'!>V't i,j,A; = 0,l,...,5, (10) 

where s is a natural parameter (time or length of the geodesic trajectory), T^- = 
\g kl (jj£ + & - §|f) is a Cristoffel symbol. Moreover x\ s = % and x\ ss = 
The system of differential equations (TTOT) is solved for the initial conditions: 

zj = x*(-oo), ±i = x\ a (-oo), (11) 



for any value of the natural parameter s from which the geodesic trajectory x l (s) 
and the geodesic velocity x l s (s) are defined. Using the relations in Eqs. fl5]) and (TTUj) 
it is not complicated to obtain the following system of equations: 




where x{ x °i xl i x2 ) = ^ 9oo(%° , x 1 , % 2 ) , the I is total angular momentum of three- 
body system. It is suitable to conduct the later calculation in the (u, v, 9) coordi- 
nates system. Because the explicit form of equation system in those coordinates is 
complicated, we don't write down them here. 

We have now formulated the reactive scattering problem in terms of classical 
dynamics on the Lagrange surface of the three-body system. Note, that the system 
has one integral of motion (overall energy E) and three degrees of freedom. 

According to Poincare (see [28}), conservative dynamical systems can have regions 
of chaotic movement in their phase space provided that they are not integrable, i.e. 
have less integrals of motion than degrees of freedom. This means that certain areas 
in phase space may show non-stability and chaos may then be observed, i.e. the 
trajectory x^(s) then becomes exponentially non-stable with respect to change of 
the initial condition x J (0): 

dx j (s)/dx j (0) ~ exp(A i s), j = 0, 1, 2, (13) 

where Xj describes the degree of instability and is called Lyapunov exponent. 
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B. Quantization of classical dynamical three-body scattering system 



Representation for regular case: 



In some coordinate systems, like the NCC system, [24j the 3-D quantum reactive 
scattering problem can be treated in the same way as an inelastic single- arrangement 
problem. The overall wavefunction of the three body system can be represented: 



J2G { ^ fg (u)E nG) (v,u)e ]KI (e), g=(n 3 K), (14) 



"J 



where (n,j,K) is a set of quantum numbers, OjK{@) is a normalized associated 
Legendre polynomial, E n (j)(v;u) is the vibrational part of the wavefunction and 
satisfies the equation: 



h 2 d 2 jjffij + i) 

2 + U (u, v) H 



2/i dv 2 



2/jv 2 



l n (j){v; u) = e n{j) (u) H n(i) 0; u), (15) 



where U(u,v) = U(u,v,6)\ 0=Q —U e ff(u,v). The function U(u, v, @)\ e=0 describes the 
potential energy of the collinear collision. U e ff(u,v) is an effective potential: 



1 / drj 



1 d 2 rj 5 / drj 



+ 



V (u,v) = [l + K{u)v]^-. (16) 



Ue ff( u ' v > Ar, 2 \dv) 2rfdu 2 ' Ar] 4 \du 
In eqn. (116j) . i^(it) is the curvature of tt and s is the length along the curve Q^-: 

[^(^)]" 3/2 



#(u) = 2a- 



ds 



{q c i-qe q f du b + a/(qi-q-) 



2 ' 



b-a/{q{- q eq )' 



Note, that the scattering function G^tH (u) satisfies the following equation: 



(17) 



d 2 

0„>ri.—r-z + 2 



Id 1 drj 



njK'g 

Id 2 2d?] d 



2/4 / 2 



n n du 2 ^ \du r] du / n , n du ^ \ du 2 rj du du / , n ^ h 2 \ 

h 2 ](j +i 



E-E 



JK 



- e n(j)( U ) + U(U,V) + 
,2 



&3'j5k>k - j^(v 2 U«,j{u, v)^ r _ 5 K , 



2{iv 2 

8 K '+1K C Jj(K-l) + s k'-ikCjj {R+1) 
where Cf jK = c^ K c^ K , and moreover :c^ K = [J(J+1) — K(K±l^ 



K 



G[% (u) = 0, 



+ ( 72 ) ,Ji'i 
I/O i n'n 



,j- K — ^J K ^j K , anumui C u» C i.u Jif — y.j T x; "iv , E JK (u,v) = 

(l/2)h 2 jj~ 1 qQ 2 (u,v)(^J(J + 1) — 2K 2 ). The summations over the repeating index n' 
and j ' are implied and we use the following notations for the matrix elements: 



(f( U ))nn> 



^n(j){v, u)f(u, v)E* n>{j) (v; u)dv, 

o 

O jK (6) U(u, v, 9) Qf K (0) sin OdB. 



(19) 
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Note, that the solution of Eq. (JT81) is must satisfy the asymptotic condition: 

lim V Gtyf K , (u) = -= exp u) <W%"W- (20) 

u^—oo 1 ' J K V 27T 

n'j ' 

The exact 5-matrix elements can be constructed in terms of stationary overall and 
asymptotic wavef unctions, considering that the variable u plays the role of a timing 
■parameter ( which later will be called internal time) (24j : 

Sn'j'K'^njK(E) = ^Prfji/Pnj J™^ X) G^tf/ e H Wfin' (w)Ajx'^jK', (21) 

nj 

where W n / B (u) = (S fi (j)(u; ^nj"^^))^ k- jK ,^ jK , = (Qj K ,(6)Q jK ,(6)) 9 = 5jj. 
The expression for the ^-matrix elements in eqn. ff2TT) can be simplified, if we take 
as basis the functions 'B n ^(v;u), which in the limit u — > +00 coincide with the 

orthonormal basis of the (out) asymptotic wavefunctions U^^(v). In this case we 
get the simplification \im u ^ +00 W n 'n(u) = 5 n 'n and the following expression holds for 
the S'-matrix elements: 

S nTK ^ njK (E) = y/pt T /p^Gy e J (E; +00), g = (njK). (22) 

Representation for chaotic case: It is well known that some chemical reactions, 
especially when highly excited, exhibit quantum chaotic behavior 0, H], i.e., the 
statistical properties of eigen-energies and eigen-vectors are very similar to those of 



random matrix systems [30|, [31 ] 



For systems which are not too quantum mechanical in nature, the quantum prob- 
ability current is localized along the classical trajectory. In the chaotic case, these 
classical trajectories diverge exponentially from each other and from the quantum 
current tubes too. This results in serious difficulties in describing chaotic reac- 
tive quantum processes in terms of standard quantum representations. In order 
to overcome the this problem, a new quantization method bases on the quasi- 
classical approach has been proposed [l9( for the three-body system. The idea is 
to carry out the quantization on separate classical trajectory tubes 3?(x 3 (s)), where 
x 3 (s) = x 3 (w(s), v, 9) and u(s) is the solution of the geodesic equations (fl2"l) . which 
varies along the curve Sfj/ and is called as a generating trajectory (recall that it 
has the meaning of internal time). Every solution u(s) generates some topological 
trajectory tube, which can be described by the Schrodinger equation, which for the 
present case means Eq. ( fl8l) . The summed contribution of all such tubes gives the 
whole quantum picture. 

The goal of the scattering problem is the calculation of the probability ampli- 
tudes for transitions between different asymptotic states. In mathematical language 
this corresponds to the determination of the total mathematical expectation of the 
elementary quantum process in the three-body system. In the classical case the 
solution u(s) depends on the initial scattering phase <p = 2ir{u/L}, as does then the 
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overall wavefunction and S- matrix elements. Here L is a some period and {•} de- 
scribes a fractional part of the function. This implies that the transition amplitude 

I 1 2 

A g r g = \S e r<- e (<p] E)\ must be averaged over the (p phase distribution: 

W e >^ e (E)- =— ; — , (26) 

J <r(ip;E)d<p 

where o~(ip; E) is the distribution of classical trajectories which will be determined 
later in the section III. In the case when the chaotic regions in phase space of classical 
system are smaller than the elementary quantum cell h (N is the dimension of 
configuration space) the transition amplitude A g i e (ip ; E) is independent from if. 



III. NUMERICAL EXPERIMENT 

Numerical calculations are here made for the collinear reaction Li + (EH) n — > 
(LiFH)* — > (LiF) m + H. The LEPS type potential energy surface of Carter and 



Murrell for this reaction was used [32[. The classical trajectory study was performed 
by solving eqn. ( |T2l) for a total angular momentum quantum number (j = 0, and 
fixing the NCC angle (Jacobi angle) to 6 = 0). In Fig. 2, three generating trajectories 
(or internal time) u(s) and their corresponding v(u) graphs are shown for different 
initial phases <p of the trajectories. It is seen that the generating trajectories behave 
quite differently depending on the initial phase ip for fixed energy E. Panel a) 
in Fig. 2 shows a direct exchange reaction to which corresponds a monotone, but 
not uniformly changing, internal time (as a function of the natural parameter s 
(usual time). Panel b) shows a non-reactive trajectory to which corresponds non- 
monotone internal time. In panel c) the geodesic trajectory again describes the 
exchange reaction which here goes via a resonance (ABC)* state and for which the 
dependence of u on the parameter s is complicated. 

Now the main task is the investigation of the behavior of the geodesic trajectory 
flow. Numerical calculations shows, that for initial values corresponding to the 
chaotic regions mentioned above, the main Lyapunov exponent is positive and grows 
fast. The last fact points to exponential divergence of geodesic trajectories. 

In Fig. 3, the white points in initial parameter space correspond to the transition 
from the reactant (-R^ n ) subspace) to product regions (-RL^) subspace), while the 
black points correspond to the reflection back to the product region. The distribu- 
tion of black and white points depend on energy and initial phase, for fixed initial 
vibrational coordinate Vq, and shows an irregular behavior. Recall that Vq is an av- 
erage equilibrium distance between bound particles B and C in the ground (n = 0, 
where n-is a vibrational quantum number) state. Note that qualitatively the same 
picture we get for v n (equilibrium distance on excited quantum stat n). One can 
see from the results of calculations that the structure of chaotic behavior region is 
self-similar with respect to scale transformation. 

Let us consider the influence of chaotic behavior of the classical problem in Fig.4, 
which show the dependence on energy E of over-barrier transition probabilities in 
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a) 



b) 
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FIG. 2: 



Geodesic trajectories and internal times depending on natural parameter 



s. 



the Li + FH system for fixed phases <fi and equilibrium distance v n . It can be seen 
that a small change in initial phase significantly changes the dependencies. In this 
connection the difficult problem arises to find a measure for the space (map) of passed 
through and reflected back geodesic trajectories. To calculate the probability for a 
specific quantum transition at an energy E, one has to average the corresponding 
quantum probability with respect to (E, ip) within the range [AE, Acp = 2ir], where 
AE is a small interval of energies near E and Atp = 2ir is the period of the values for 
the initial vibrational phase. The procedure of the averaging consists of that square 
[AE x 2tt} is divided by N ^> 1 rectangles, each of them having some phase point 
(fi inside. Then each rectangle is subdivided by the grid with Mj = [Zj x kj\ nodes, U 
and ki being the number of breaking points for AE and 2ir/N intervals respectively. 
Probability for geodesic trajectory {generating trajectory) to pass through the i-th 
rectangle is calculated by the formula: 



a{cpi, E) = lim Nij Mi 



(24) 



where Ni counts how many times the generating trajectory passes through into 
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FIG. 4: (a) - dependencies of transition probabilities Aoo, Aoi and A02 on energy E for 
fixed phase (p; (b) - the same dependencies, but calculated for the other (slightly differing 
from the first one) fixed phase (p a — fb = 10 -5 . 

subspace Rf in \ ■ Exchange reaction probability is then calculated as a limit of sum: 
w y I T,?=i cr (<Pi)\ S n^m(<Pi,E)\ 2 \ = fy n cr(<p)A nm (<p,E)d<p 



Particularly using (|25|) for reacting system Li + (FH) — > (LiFH)* — > (LiF) +H we 
can calculate transition probabilities see Fig.5. 
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FIG. 5: The transition probabilities after averaging with respect on phase distribution o~((p). 



IV. CONCLUSION 

The 3D quantum theory Eq.s f|T2|) and f|T8|) was constructed, which described 
the classical permissible P 2 (u,v,8) > reactive scattering in the three-body sys- 
tem with taking into account the influence of classical non-integrability on the 
quantum dynamics. By means of numerical modelling of collinear reacting sys- 
tem Li + (FH) n — > (LiFH)* — > (LiF) m + H it was shown that when chaotic region 
in phase space of classical system is larger than the quantum sell volume in the 
corresponding quantum system chaos is generated too. The calculations which are 
made for reacting collinear systems O + N 2 , N + N 2 , N + 2 , N + 2 show that 
the classical chaos is not enough developed and can't generate the quantum chaos. 
However not excepting that these systems can be quantum chaotic in the 3D scat- 
tering case. Moreover we are sure that all 3D three-body scattering systems in more 
or less degrees are chaotic. 

Finally it is necessary to note that developed approach to give a chance by defined 
conditions pass to R, P and Q regions of motion. When classical chaos is absent 
or developed insufficiently strong this approach coincides with standard quantum 
representation. 
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